Impurity transport in ITG and TE mode dominated turbulence 

A. Skyman 1 , H. Nordman 1 , P. Strand 1 , F. Jenko 2 , F. Merz 2 

1 Euratom-VR Association, Department of Radio and Space Science, Chalmers University of 

Technology, SE-412 96 Goteborg, Sweden. 
2 Max-Planck-Institut fiir Plasmaphysik EURATOM-IPP, D-85748 Garching, Germany. 



O 

(N 

O 

<D 

Q 

o 

(N 

s 

O 

>V 
^ . 
l Qh. 

> 

oo 

o 

m 

o 

> 

IS 





























Slope /), ; 

\ 







































Introduction 

The transport properties of impurities is of high relevance 
for the performance and optimisation of magnetic fu- 
sion devices. For instance, if impurities from the plasma- £ 
facing surfaces accumulate in the core, wall-impurities of 
relatively low density suffice to dilute the plasma and lead 
to unacceptable energy losses in the form of radiation. 

In the present study, turbulent impurity transport in Figure 1: musfratio n of PF and the lin- 
Deuterium tokamak plasmas, driven by Ion Temperature earity of T z (Vn z ); ITG dominated quasi- 
Gradient (ITG) and Trapped Electron (TE) modes, has linear GENE result for Ne with k P = 3 
been investigated using fluid and gyrokinetic models. The impurity diffusivity (D z ) and con- 
vective velocity (Vz) are calculated, and from these the zero-flux peaking factor (PF ) is de- 
rived. This quantity expresses the impurity density gradient at which the convective and diffu- 
sive transport of impurities are exactly balanced. The sign of PF is of special interest, as it 
determines whether the impurities are subject to an inward (PF > 0) or outward (PF < 0) 
pinch. 

Quasilinear results obtained from the GENE code flU 13 are compared with two-fluid results 
for both ITG and TE mode dominated turbulence. Scalings of PF with impurity charge (Z) 
and various plasma parameters, such as magnetic shear (s), are studied. Of particular interest 
are conditions favouring an outward convective impurity flux. 

Theoretical background 

The transport of a trace impurity species can locally be described by a diffusive and a convective 
part. The former is characterized by the diffusion coefficient D z , the latter by a convective 
velocity or "pinch" V z , see equation (Q~|) 01. From these, the zero flux peaking factor is defined 
as PF = — ^p-| r=0 , see figure[TJ PF is important in reactor design, as it quantifies the 
balance of convective and diffusive transport. This can be seen from equation (Q~|), where T z is 
the impurity flux, n z the density of the impurity species and R the major radius of the tokamak. 



For the domain studied - a narrow flux tube - the gradient of the impurity density is constant: 
Vnz/nz = 1/L nz . Setting T z = in equation (OQ) yields the interpretation of PF as the 
gradient of zero impurity flux. 



.Rr r 

Tz = -DzVnz + n z Vz^ = -D z — + RV Z (1) 

n z L nz 

Fluid model 

Though the main results presented in this study have been obtained using quasilinear gyrokinetic 
simulations, their physical meaning is interpreted by comparing with the Weiland multi-fluid 
model 0. The fluid equations for each included species (j = i, te, Z , representing Deuterium 
ions, trapped electrons, and trace impurities) are: 

OTL ■ 

-^ + V-( W ) = (2) 

dv\\i z 

mi,zni,z dt + V|| {n it zT it z) + n it zeV\\(p = (3) 



3 d^Z" 1 

^ ] ^ + n J T ] W-v ] + y-q J = Q (4) 



Here q^ is the diamagnetic heat flux, and Vj is the sum of the E x B, diamagnetic drift, po- 
larization drift, and stress-tensor drift velocities. To solve the equations, it is assumed that is 
the only heat flux for all species, that passing electrons are adiabatic, and that quasineutraility 
(equation ©) holds. Going to the trace limit for the impurities, i.e. letting Zf z — >■ in equa- 
tion ©, an eigenvalue equation for ITG and TE modes is obtained. The impurity particle flux 
in equation (QQ) is then obtained from T n j = (5rijt; BxB ), where the averaging is performed over 
all unstable modes for a fixed length scale kp of the turbulence. 

^ = (1 - Z h ) ^ + Zh^-, h = - (5) 
n e m n z n e 

Quasilinear gyrokinetic simulations 

GENE is a parallel gyrokinetic code employing a fixed grid in five dimensional phase space and 

a flux-tube geometry [1J. The simulations were performed on the HPC-FF clusteiQ with GENE 

running in eigenvalue mode. Growth rates and impurity fluxes were thus computed for ITG and 

TE mode dominated cases, for which a number of parameters were varied and trends observed. 

The main parameters used are presented in table [TJ 

'HPC-FF (High Performance Computing For Fusion) is an EFDA funded computer situated at Forschungs- 
zentrum Jiilich. Germany, dedicated to fusion research 



Table 1 : Parameters used in all simulations 





ITG: 


TEM: 


T D /T e : 


1.0 


1.0 


s: 


0.8 


0.8 


<Jo: 


1.4 


1.4 


e: 


0.14 


0.14 


R/ l t d , R/ Lt z ■ 


7.0 


3.0 


R/Ltc ■ 


3.0 


7.0 


N x x N ky x N z : 


5 x 1 x 24 


4 x 1 x 24 


N V{] x iV M : 


64 x 12 


64 x 12 



Results 

Impurity charge Z: The main results obtained are the scalings of the peaking factor with the 
charge of the impurity species. These are presented in figures [2(a)] and |2(b~)| showing ITG and 
TE mode dominated turbulence respectively. 
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(a) ITG mode dominated case (b) TEM dominated case 

Figure 2: Scalings of PFq with impurity charge Z; quasilinear GENE and fluid results 



The difference between figure [2(a)] and [2(b)] can be understood from the properties of the con- 

1 11 and a parallel impurity 



vective velocity in (Q])- Vz contains a fhermodiffusive term Vr z 



compression term V p 



-/,- 



r . The former is generally outward (Vt z > 0) for ITG 



Vz A z n \\ A z qf 

and inward (Vt z < 0) for TE mode dominated transport, whereas for the latter the opposite is 
generally the case. 



Magnetic shear s: The effect of magnetic shear on the peaking factor is shown in figures |3(a) 
and |3(b)[ It is worth noting that a flux reversal, i.e. a change of sign in PF , owing to a change 
in sign of Vz, occurs for negative s for Z > 6 in the TE mode dominated case, indicating a net 
outward transport of the heavier elements. Similar trends are not seen in fluid simulations, and 
this warrants further investigation. 



(a) ITG mode dominated case (b) TEM dominated case 

Figure 3: Scalings of PFq with magnetic shearing s; quasilinear GENE results with kp = 0.3 



Other parameters: Scans of the dependence PF on other parameters, such as kp and Lt, 
have also been carried out. The results are similar to those reported in Q, O and Q respec- 
tively. In most cases, only a weak dependence of PF is observed. 

Conclusions and outlook 

Quasilinear GENE simulations and fluid results show that peaking factor increases with im- 
purity charge Z for ITG mode dominated transport, whereas the opposite holds for TE mode 
dominated transport. In both cases PF saturates for high Z. 

For magnetic shear, a flux reversal is observed for negative magnetic shear in the TEM dom- 
inated case. This is not seen in fluid simulations, and will be a focus of future studies. 

For other parameters investigated, weak scalings for PF are observed, in agreement with 
previous work. 
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